Repulsion-to-attraction transition in correlated electron systems triggered by a monocycle pulse 
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We study the time evolution of the Hubbard model driven by a half-cycle or monocycle pulsed electric field 
F(t) using the nonequilibrium dynamical mean-field theory. We find that for properly chosen pulse shapes the 
electron-electron interaction can be effectively and permanently switched from repulsive to attractive if there is 
no energy dissipation. The physics behind the interaction conversion is a nonadiabatic shift S of the population 
in momentum space. When 6 ^ n, the shifted population relaxes to a negative-temperature state, which leads to 
the interaction switching. Due to electron correlation effects S deviates from the dynamical phase </> = J dtF(t), 
which enables the seemingly counterintuitive repulsion-to-attraction transition by a monocycle pulse with </> = 0. 

PACS numbers: 71.10.Fd, 03.65.Vf, 03.75.Ss 



I. INTRODUCTION 

Controlling the interparticle interactions in a correlated 
electron system by applying intense laser fields is a challeng- 
ing and exciting perspective, which may lead to states of mat- 
ter that do not exist in equilibrium. For example, if one could 
effectively change the electron-electron interaction from the 
original Coulomb repulsion to an attraction, this may induce 
an i-wave superconducting state with very high transition 
temperature (optimally « 0.1 bandwidth),''^ or the BCS-BEC 
crossoveri^il It will also enable us to study phenomena char- 
acteristic of nonequilibrium quantum systems, such as tran- 
sient states after an interaction quench."*"^ The control of inter- 
particle interactions is in fact possible in cold-atom systems,^ 
where one can manipulate the interaction in a wide range from 
repulsive to attractive using the Feshbach resonance that dom- 
inates the scattering length,^ but such a technique cannot be 
applied to electron systems. 

One way to control the interaction is to create a popula- 
tion inversion in metallic bands corresponding to a negative- 
temperature (T) state. This implies an effective switching 
of the interaction from repulsive to attractive, since a density 
matrix e^^^^ for a Hamiltonian H with temperature T < Q 
corresponds to the one for the inverted H with —T> 0."^- 
While a (partial) population inversion itself is a common phe- 
nomenon (e.g., in laser productions), the interaction conver- 
sion is a genuine coiTelation effect in nonequilibrium. Ideally, 
the laser fields that drive the system should be "pulsed" waves 
since, first, the available intensity is generally much higher for 
ultrafast pulses'^"'^ than for continuous-wave lasers, and sec- 
ond, continued heating can be avoided. These considerations 
raise a fundamental question: can irradiation by a single-cycle 
pulse put a system into a negative-T state that survives for a 
long time after the pulse? 

In this paper, we show that it is possible to induce a pop- 
ulation inversion in metallic systems using a properly shaped 
monocycle or half-cycle pulse, and that in the absence of en- 
ergy dissipation, the system will thermalize in the negative- 
T state after the pulse. By solving the driven Hubbard 
model with the nonequilibrium dynamical mean-field theory 
(DMFT)r^iii we will demonstrate that pulse fields F{t) with 



proper asymmetry between the positive [F(t) > 0] and neg- 
ative [F(t) < 0] parts trigger a repulsion-to-attraction transi- 
tion. Such asymmetric pulses can readily be generated thanks 
to the recent progress in laser techniques, '^^"^ while their po- 
tential application to correlated systems has remained unex- 
plored, in contrast to symmetric, many-cycle pnlsesM^ 

Our strategy is to induce a nonadiabatic shift (denoted 
by 6) in the momentum distribution of the electrons by the 
asymmetric monocycle pulse (see Fig.[T]|. If we can achieve 
5 ^ n (half of the Brillouin zone), the system is brought to a 
negative-T state [Fig.fTfc)], which amounts to a change of the 
interaction from repulsive to attractive. In a one-body picture, 
one expects that each electron acquires from the pulse field a 
dynamical phase 



- / dtF{t) 



(1) 



with e the elementary charge and a the lattice constant (here- 
after we set e = fl = ft = 1). This causes a momentum shift 
A: — > k + (f), so that we simply have 5 - 4'- An immediate 
question is: can a monocycle pulse with j dtF{t) - (as dic- 
tated by Maxwell's equation^-*) induce a nontrivial shift of the 
population? We show that it is in fact possible in interacting 
systems, where the nonadiabatic shift 6 exhibits a clear de- 
viation from 4> due to coiTelation effects. This allows us to 
achieve 5 ^ n even when J dtF{t) oc = by choosing the 
pulse shape appropriately. We reveal conditions for the pulse 
shape that lead to the population inversion, and construct a 
'phase diagram' for the pulse-driven Hubbard model. We em- 
phasize that the interaction conversion robustly persists after 




FIG. 1. (Color online) Schematic band pictures for the pulse-induced 
phase shift, (a) Initial system in equilibrium, (b) right after the pulse 
excitation, with the population shifted in momentum space by S, and 
(c) the system finally thermalized with a negative temperature when 
(5~7r. 
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the pulse has passed, at least in an isolated system without 
energy dissipation. This contrasts with the previously pro- 
posed scenario for the repulsion-to-attraction transition using 
continuous-wave fields»i^ 



II. MODEL AND METHOD 

We take, as the simplest model for correlated electrons, the 
single-band Hubbard model driven by an electric field with 
the Hamiltonian 



(a) 



(b) 



t 



FIG. 2. (Color online) Schematic temporal profiles of a half-cycle 
pulse (a) and a monocycle pulse (b) for which a controls the asym- 
metry. 



exp 



i dr ■ A{t) 



(2) 



where tij is the hopping between sites at and Rj, the electric 
field F{t) - —dA(t)/dt is expressed in terms of the vector 
potential A(t), and (c) creates (annihilates) an electron. For 
the interaction we take the particle-hole symmetric form 



(3) 



where U(> 0) is the repulsive Coulomb interaction with 
n/o- - clj-Cia-. We apply a pulsed wave at f = 0, and switch 
off the field at f = t. For the DMFT, we consider a hy- 
percubic lattice with the Gaussian density of states D(e) - 
^ and apply the field in the diagonal direction 
with F{t) - F{t){l, 1, . . .). The band is assumed to be half- 
filled. Throughout the paper, we use the bandwidth W as the 
unit of energy, and take the initial temperature to be T = 0.1. 



III. RESULTS 



A. Noninteracting system 



Let us start with the noninteracting system. We fo- 
cus on the momentum distribution defined by f(k, t) - 
-iGt(t,t) = -/G<^^(,)(f,f), where G^itj') [G^(f,f')] is the 
(gauge-invariant^"*) lesser Green function. For the noninter- 
acting system, the lesser Green function is given by 



Gokit, t') = ifoiek) exp ( / dt£t-A(r) 



(4) 



where /o(e) = 1/(6*^^^ -H 1) is the Fermi distribution, and 
eic the band dispersion. After the pulse excitation (f > t), 
the momentum distribution becomes f(k, t) - fo(ek-^) with 
^ - -A(t) - 0(1, 1, . . .). Note that the effect of the pulse 
field on the final state amounts to a momentum shift ([T]). 
For a 7T shift ((/> ~ n), the electrons occupy the band top with 
f{k, t) ~ /()(— ft), which is characterized by an effective tem- 
perature Teff = — r < 0. Thus a n shift (which may be viewed 
as a partial Bloch oscillation) is the condition that leads to a 
negative-T state in the noninteracting system. 



B. Interacting system 

Now let us move on to the interacting case. There, we can 
identify the repulsion-to-attraction transition from the total 
energy Etotif) - {H{t)): After the pulse excitation, a (nonin- 
tegrable) isolated system is supposed to approach a thermal- 
ized state^^ with some effective temperature T^a and a total 
energy Eiotij) (which is conserved after the pulse is over at 
t - t). a thermal state with a positive temperature always 
gives £(01 < at half filling for the interaction term ([3]), while 
one with a negative temperature gives fitot > 0. This suggests 
that the total energy plays the role of an "order parameter" for 
the repulsion-to-attraction transition. If and only if Eiotij) > 
the system arrives at a negative-T state {T^s < 0), in which 
case the density matrix is given by 



p cx exp 



exp 



1 
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(5) 



Here we have introduced a gauge transformation 

,-<(0+;r)(I,I„ 



Ci,j 



with Hi„t(-U) = -t/E/ («'t - 5) («4 - 5) (with hi^ = 
cj^c/o-) to cancel the phase shift in the kinetic energy. The 
above equation implies that the state can be viewed as a 
thermal state with a positive T and an attractive interaction 
~U < Q. This is the basic mechanism behind the repulsion- 
to-attraction transition driven by the pulse. Note that the con- 
dition for the interacting system is Teff < 0, as opposed to the 
noninteracting counterpart ~ tt. 

To make our argument more precise, we consider two types 
of pulses. One is a half-cycle pulse [Fig.|2ja)], and the other 
is a monocycle pulse [Fig.|2tb)]: 



^'half-cycle(f) - 



monocycle 
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Here A controls the amplitude of the pulse, s{x) (> 0) is a pulse 
shape function that has support in < jc < 1 with Jy dxs(x) - 
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FIG. 3. (Color online) Phase diagram in the (U,A) plane for the 
Hubbard model driven by the half-cycle cosine pulse. The regions 
surrounded by solid curves represent the pulse-induced, attractively 
interacting phase, while the dashed curves are loci of the shift S = 
(2n + l)n. The inset shows the phase diagram in the (U, t) plane for 

A = 7T. 



I, and a (0 < Q- < 1) controls the asymmetry of the monocycle 

letric case [F{t — 
A (0 = 0) for the 



pulse, with a - ^ corresponding to the symmetric case [F{t 



t) - —F{ty\. The dynamical phase ([TJ is 
half- (mono-)cycle pulse. 
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C. Half -cycle pulse 

We first consider the half-cycle pulse (|6]l. The simplest case 
is the limit t — ^ 0, corresponding to a delta-function pulse 
[F{t) — > A6{t)]. In this case, the momentum shift is 5 = 0, so 
that the order parameter reads 

Etotir) = £'kin(0) cos (/) - iJ(0) sin + £'int(0), 

where E]^in(t), J{t), Eint{t) are the kinetic energy, current, and 
interaction energy at time f, respectively. Since 7(0) = in 
the initial state and <p = A for the half-cycle pulse (|6]l, we 
can identify the condition for the repulsion-to-attraction tran- 
sition, 

£kin(0)cosA + £i„,(0)>0, (8) 

which is completely determined by the equilibrium state. The 
criterion is quite general, so should be applicable to systems 
with any fillings in any dimensions if one puts the origin of 
the energy to be the one atT - ±oo. In Fig. |3] we show the 
attractively interacting regions by the hashed areas. Attractive 
regions appear periodically in the amplitude A of the pulse as 
a series of lobes around A ~ (2« -i- l);r (« = 0, 1, 2, . . .). Each 



FIG. 4. (Color online) Time evolution of the double occupancy (a) 
and the current (b) after the delta-function-pulse excitation, which 
are compared with the interaction quench (U —U) with U = 
1 . Each arrow indicates the double occupancy in the corresponding 
thermal state with the same total energy. 



lobe has the tip at Uc — 2.186, which turns out to be smaller 
than the critical U for the Mott transition,'^ so that the tran- 
sition always occurs in the metallic regime. The repulsion- 
to-attraction conversion is obviously distinct from a heating 
effect, since it appears and disappears repeatedly as one in- 
creases the amplitude A of the pulse field. 

To study how the system evolves in time, we have numeri- 
cally solved the model (|2| with the nonequilibrium DMFT. — 
As an impurity solver for DMFT, we mainly employ the 
continuous-time quantum Monte Carlo method^^ with the 
weak-coupling expansion generalized to nonequilibrium.^^ To 
capture the long-time (f > 20) behavior in the very weak- 
coupling regime (U < 1.2) we use the iterative perturbation 
theory, "'•^'^ which is known to give quite accurate results up to 
a long time^^ for such small U at half filling. 

In Fig. Ufa) we show how the double occupancy d(t) = 
{n^n^{t)), a measure of the interaction, evolves after the delta- 
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(a) A = 0.8;r 



1 




(b) A = ;r 



1 




FIG. 5. (Color online) Snapshots of the corresponding momen- 
tum distribution f{k, (black curves) at ? = 0, 0.8, 1.6, . . . , 16 are 
shown for A = 0.8;r (a) and A = ;r (b). They converge to thermal dis- 
tributions with Teff = —0.725 and Teff = —0.446 (thick red curves), 
respectively. Dashed curves represent the initial distributions. The 
horizontal arrows indicate the pulse-induced phase shift 6 = (f> = A. 



function pulse in an initially repulsive system {U — 1). We 
notice that for A > O.Stt d(t) shoots well beyond the nonin- 
teracting value d = {n^){ni) - 0.25, which implies that the 
electrons do indeed start to attract each other after the pulse, 
as predicted from the criterion The repulsion-to-attraction 
transition is "perfect" for A = tt, where the temporal evolu- 
tion of d is found to accurately agree with that for the inter- 
action quench, U ^ —U [Fig.|3a)]. For this pulse", the 
shift of the momentum just changes the sign of the hopping 
(ek ek+if, = —ek), which is known to be equivalent to inter- 
action quenchji^ 

Remarkably, after the pulse excitation d(t) relaxes quickly 
to a steady state (f < 3). We have confirmed that it converges 
to the thermal value t/th [indicated by arrows in Fig. Ufa)] for 
the equilibrium state having the same istot- For A > O.Stt the 
coiTesponding temperature (Teff) of the thermal state is nega- 
tive since Etot > 0. Note that dth is a nonmonotonic function 
of temperature, so that d{t) decreases in time for A - Q.2n. 

On the other hand, the current 7(f) [Fig.|3b)] generated by 
the momentum shift k k + (p decays more slowly (t < 15) 
for < A < TT than d{t). The slow relaxation is also seen in 
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FIG. 6. (Color online) The total energies at f = r for the system with 
U = \ driven by monocycle pulses with A = tx, t = plotted for 
various pulse shapes as a function of the asymmetry parameter cy. 

the momentum distributions f{k, [Fig. 13 a) for A - O.Stt 
and (b) for A - n]. A similar separation of the relaxation 
times of d( f) and f{k, t) has been observed in the interaction 
quench,*" and was attributed to the existence of a "prethermal- 
ized" state. ^ Here the slow decay becomes particularly evident 
when the momentum shift is not perfect (i.e., A 4^ n). In this 
case the system needs to adjust the momentum shift to tt to 
achieve a thermal state. Since the relaxation involves a viola- 
tion of the momentum conservation by Umklapp scattering, it 
takes longer than the relaxation of the double occupancy via 
particle-hole annihilations. The distributions eventually relax 
to thermal states with Tefj < [red curves in Fig.jSfa), (b)]. 

So far we have examined the delta-function pulse (r — ^ 0). 
To be more realistic it is important to evaluate the effect of t 
on the transition. Here we take, as an example, a half-cycle 
pulse with 

s{x) - \ — cos(2ttx), (9) 

which we call the "cosine pulse". We can again use £tot(T) 
as an indicator for the interaction conversion. In the inset 
of Fig. [3j we show how the critical interaction iUc) of the 
repulsion-to-attraction transition induced by the cosine pulse 
with A - n depends on r. For r < 10, Uc{t) rapidly falls 
off from Ucij - 0), while for larger t it decays to zero very 
slowly. In the adiabatic limit (t — > oo) Eiotij) £101(0) < 
for U > and the repulsion-to-attraction transition naturally 
disappears. 

For general A and non-zero t of the cosine pulse, Fig.[3]de- 
picts the phase diagram for the pulse-driven Hubbard model. 
The attractive regions now deform in a characteristic manner, 
i.e., the tips of the lobes bend toward larger A, which becomes 
drastic for r - 40. The deformation of the phase diagram im- 
plies a rather counterintuitive fact: for t = 40 and U ^ 0.9 
the repulsion-to-attraction conversion occurs even for a triv- 
ial phase (p - A = 2tt. This suggests that the effective phase 
shift 6 that the correlated system acquires is not equal to <p 
for nonzero t. Using the extremal points of the total energy 
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satisfying 5£'tot(T")/5A = (dashed curves in Fig.|3]l as an es- 
timate for 6 = (2n + l)7r, we can see the large deviations of 
6 from = A = (2« + l);r as U and t grow. We attribute 
this to a correlation effect: During irradiation with the pulse 
the electrons scatter with each other, which causes broadening 
of the momentum distribution. Consequently the shift in the 
momentum is suppressed, and 6 becomes smaller than (p. 



case of U - I and r = 40, the larger the peak amplitude 
[i (i) - |,2, for the cosine^, cosine, sine, semicir- 

cular, rectangular pulse, respectively], the larger the £'tot('r) 
around a - ^. For even shaper pulse shapes (cosine^*, . . .), 
Etot(T) at Q- = ^ starts to decrease, so that the cosine^ pulse is 
an optimal shape in this case. 



D. Monocycle pulse 

The result that 5 4' suggests an experimentally much 
simpler way to induce the repulsion-to-attraction transition 
{5 ~ n) by monocycle pulses (|7]i with = 0. The basic idea is 
the following: since the effect of a half-cycle pulse very much 
depends on its width, we can suitably choose the widths of the 
first and second half cycles of a monocycle pulse so that the 
total phase shift is c± n. In Fig.|6l we plot the total energy at 
t - T for various types of pulse shapes. The shape function of 
each pulse is defined by 

— cos(27rx)]^ cosine^ pulse, 

I sminx) sine pulse, 

^ -\/l — (2x — 1)2 semicircular pulse, 

1 rectangular pulse, 

for < X < 1 . six) of the cosine pulse is defined as before 
[Eq. (|9j]. From the above argument, the asymmetry, here rep- 
resented by a, should be important, and we can indeed see that 
Etoti^) becomes positive (implying a repulsion-to-attraction 
transition) as soon as we go sufliciently away from the sym- 
metric pulse form {a - j). To be more precise, the momen- 
tum shift induced by the first half cycle does not cancel the 
one induced by the second half cycle when the mono-pulse 
is sufliciently asymmetric. Note that the cosine^ pulse gives 
£tot(T) > even at a = |, although its value is very small. 
There is a slight difference in £'tot(T) for the cases with a < ^ 
(i.e., the sharp pulse comes first, followed by the broad one) 
and a > ^, implying that the former is more suitable than 
the latter to induce the attractive interaction. The order pa- 
rameter £'tot(T) also depends on the shape of the pulse. In the 



IV. DISCUSSION AND OUTLOOK 

Finally, let us discuss the experimental feasibility of 
the pulse-induced repulsion-to-attraction transition proposed 
here. Asymmetric monocycle pulses with qualitative features 
comparable to the shapes considered here can be generated 
experimentally.'" One way to detect the negative-T state is to 
measure the time-resolved dc or optical conductivity, which 
will become negative after the pulse irradiation due to energy 
gain. Another possibility is to measure momentum-resolved 
photoemission spectra, which can detect the shift in the mo- 
mentum distribution [Fig.|5la), (b)]. We require the time reso- 
lution of the measurement to be fine enough that it can detect 
the population-inverted state before it relaxes to a more sta- 
ble state through energy dissipation. The dissipation typically 
occurs due to phonons whose time scale is of the order of 
0.1 — 1 ps,''^ which allows one to access the negative-T state 
using current ultrafast laser techniques with a resolution ~ 10 
fs.''* Materials that have a metallic band at the Fermi energy, 
well separated from the other bands, are suitable candidates 
because of the absence of interband transitions that destabi- 
lize the population-inverted state. As an example, transparent 
conductors (e.g., Sn-doped In203^^) and alkali-metal-loaded 
zeolites^"* such as sodalite^ are materials that exhibit this kind 
of band structure. 



V. ACKNOWLEDGMENTS 

We thank C. Hauri, C. Mudry and R. Shimano for stim- 
ulating discussions and S. Biermann for suggesting Sn-doped 
In203 as a potentially interesting material. This work was sup- 
ported by SNF Grant No. PP0022-118866, and by a Grant- 
in-Aid for Scientific Research on Innovative Areas "Optical 
Science of Dynamically Correlated Electrons". 



' R. Micnas, J. Ranninger, and S. Robaszkiewicz, Rev. Mod. Phys. 
62, 113 (1990). 

^ M. Keller, W. Metzner, and U. SchoUwock, Phys. Rev. Lett. 86, 
4612 (2001). 

' P. Nozieres and S. Schmitt-Rink, J. Low Temp. Phys. 59, 195 
(1985). 

* S. R. Manmana, S. Wessel, R. M. Noack, and A. Muramatsu, 

Phys. Rev. Lett. 98, 210405 (2007). 
^ M. Moeckel and S. Kehrein, Phys. Rev. Lett. 100, 175702 (2008). 

M. Eckstein, M. KoUar, and P. Werner, Phys. Rev Lett. 103, 

056403 (2009). 



^ I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 80, 885 
(2008). 

^ C. Chin, R. Grimm, P. Julienne, and E. Tiesinga, Rev. Mod. Phys. 

82, 1225 (2010). 
' E. M. Purcell and R. V. Pound, Phys. Rev. 81, 279 (1951). 

N. R Ramsey, Phys. Rev. 103, 20 (1956). 
" A. Rapp, S. Mandt, and A. Rosch, Phys. Rev. Lett. 105, 220405 

(2010). 

N. Tsuji, T. Oka, P Werner, and H. Aoki, Phys. Rev. Lett. 106, 
236401 (2011). 

" A. L. Cavalieri, N. Muller, T. Uphues, V. S. Yakovlev, A. Bal- 
tuska, B. Horvath, B. Schmidt, L. Blumel, R. Holzwarth, S. Hen- 



6 



del, et al, Nature (London) 449, 1029 (2007). 
S. Wall, D. Brida, S. R. Clark, H. P. Ehrke, D. Jaksch, A. Ardavan, 
S. Bonora, H. Uemura, Y. Takahashi, T. Hasegawa, et al., Nat. 
Phys. 7, 114 (2011). 

R. Ulbricht, E. Hendry, J. Shan, T. F. Heinz, and M. Bonn, Rev. 
Mod. Phys. 83, 543 (2011). 

A. Georges, G. Kotliar, W. Rrauth, and M. J. Rozenberg, Rev. 
Mod. Phys. 68, 13 (1996). 

" J. K. Freericks, V. M. Turkowski, and V. Zlatic, Phys. Rev. Lett. 
97, 266408 (2006). 

R. R. Jones, D. You, and P H. Bucksbaum, Phys. Rev. Lett. 70, 
1236 (1993). 

" J. Hebling, K.-L. Yeh, M. C. Hoffmann, B. Bartal, and K. A. Nel- 
son, J. Opt. Soc. Am. B 25, B6 (2008). 

^° M. Jewariya, M. Nagai, and K. Tanaka, J. Opt. Soc. Am. B 26, 
AlOl (2009). 

B. Moritz, T. P. Devereaux, and J. K. Freericks, Phys. Rev. B 81, 
165112 (2010). 

M. Eckstein and P Werner, Phys. Rev. B 84, 035122 (201 1). 
It is nevertheless possible to produce pulses containing only a 
half-cycle oscillation that gives nonzero (j). See, e.g., Ref.J^ 
^* N. Tsuji, T. Oka, and H. Aoki, Phys. Rev. B 78, 235124 (2008). 



See, e.g., M. Rigol, V. Dunjko, and M. Olshanii, Nature (London) 
452, 854 (2008), and references therein. 

E. Gull, A. J. Millis, A. L Lichtenstein, A. N. Rubtsov, M. Troyer, 
and P Werner, Rev. Mod. Phys. 83, 349 (201 1). 
" P Werner, T. Oka, and A. J. Millis, Phys. Rev. B 79, 035320 
(2009). 

A. Georges and G. Kotliar, Phys. Rev. B 45, 6479 (1992). 

M. Eckstein and P Werner, Phys. Rev. Lett. 107, 186406 (2011). 

k is defined such that ek = — cos k and = J^, ^ = sin 

(-;r < yt < n). 

C. Hauri, private communication; S. Watanabe, N. Minami, and 
R. Shimano, Opt. Exp. 19, 1528 (2011). 

L. Perfetti, P. A. Loukakos, M. Lisowski, U. Bovensiepen, 
H. Eisaki, and M. Wolf, Phys. Rev. Lett. 99, 197001 (2007). 
" O. N. Mryasov and A. J. Freeman, Phys. Rev. B 64, 233111 
(2001). 

R. Arita, T. Miyake, T. Kotani, M. v. Schilfgaarde, T. Oka, 
K. Kuroki, Y. Nozue, and H. Aoki, Phys. Rev. B 69, 195106 
(2004). 

K. Nakamura, T. Koretsune, and R. Arita, Phys. Rev. B 80, 
174420 (2009). 



